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Preface 

This  thesis  presents  my  attempts  to  model  the  complex  air  refuel- 
ing task  In  a relatively  simplistic  manner  in  order  to  obtain  some 
"first  cut"  predictions  as  to  Its  vulnerability  to  the  adverse  effects 
of  air  turbulence  and  pilot  irradiation.  The  study  was  done  at  the 
United  States  Air  Force  School  of  Aerospace  Medicine,  at  Brooks  Air 
Force  Base,  Texas,  where  a large  data  base  had  already  been  established 
for  the  work.  Several  very  broad  assumptions  and  simplifications 
were  made  to  reduce  this  problem  to  manageable  proportions  considering 
the  time  and  equipment  constraints  under  which  I worked.  As  a result, 

1 believe  that  the  trends  exhibited  by  the  data  and  its  sensitivity 
to  changes  in  the  simulation  variables  are  more  significant  than  the 
specific  times  computed  for  the  refueling  task. 

I have  tried  to  simulate  the  refueling  task  as  realistically  as 
time  and  equipment  permitted.  To  this  end  I observed  the  "flying  abil- 
ities" of  the  pilot  models  to  assure  that  they  were  not  controlling 

in  a superhtiman  fashion.  There  is  considerable  doubt  in  some  areas  as 

/ 

to  the  limitations  of  the  human  pilot.  For  example,  it  seems  uncertain 
as  to  how  rapidly  and  over  how  large  a range  the  pilot  can  be  reasonably 
ekpected  to  move  the  throttles  to  maintain  position  in  the  refueling 
task.  I tried  also  to  keep  the  execution  of  the  refueling  task  itself 
within  established  limits.  A good  example  of  this  action  is  the  setting 
of  realistic  limits  on  the  closure  rate  during  hookup. 

I wish  to  thank  the  staff  of  the  Aerodynamics  Branch  of  Boeing/ 
Wichita  for  the  hours  of  telephone  time  as  well  as  the  copies  of  tech- 
nical reports  I needed  for  this  work;  the  personnel  of  the  Biometrics 
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Section  of  the  USAF  School  of  Aerospace  Medicine  for  their  always 
eager  support  and  assistance;  Dr.  Charles  P.  Hatsell  for  his  help  In 
enabling  me  to  gain  a background  In  the  theory  behind  this  work,  and 
for  his  very  flexible  consultation  hours  during  the  period  spent  In 
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ABSTRACT 

A hypothetical  air  refueling  operation  Involving  a B-52H  bomber 
and  a KC-135  tanker  Is  modeled  on  an  analog  computer  for  the  purpose 
of  determining  Its  vulnerability  to  air  turbulence  and  pilot  gain 
decrement  due  to  Irradiation.  The  model  Includes  the  following 
elements:  (1)  the  longitudinal  small  perturbation  equations  of  motion 

for  the  B-52H  In  the  downwash  flow  field  of  the  KC-135;  (2)  pilot 
models  to  operate  the  elevator  and  throttle  of  the  B-52H;  and  (3)  a 
stochastic  gust  model  which  provides  vertical  gust  disturbances  to 
the  equations  of  motion. 

The  controlled  elements  with  and  without  downwash  are  derived 
for  the  refueling  task.  Comparison  of  these  transfer  functions  shows 
that  the  downwash  causes  Instability.  Stability  analyses  are  con- 
ducted on  the  pilot-aircraft  open-loop  transfer  functions  over  a 
range  of  pilot  rate  to  displacement  weighting,  Tl,  and  pilot  gains 
are  found  which  will  stabilize  the  system  for  each  value  of  Tj^. 

By  means  of  no-gust  hookup  simulations  the  range  of  Tl  Is  nar- 
rowed by  eliminating  those  values  of  Tl  which  yield  unrealistically 
long  hookup  times  or  poor  aircraft  control  during  hookup. 

Monte  Carlo  simulations  of  the  refueling  task  are  run  at  varying 
combinations  of  Tl,  rms  gust  Intensity,  pilot  thrust  authority,  and 
pilot  gain  decrement  due  to  Irradiation.  Trends  of  mean  refueling 
time  are  predicted.  Using  mean  refueling  time  as  a basis  for  com- 
parison, optimum  values  of  Tl  are  found  for  the  Irradiated  and  the 
non-lrradlated  pilot. 

Results  also  are  presented  on  a small  group  of  Monte  Carlo  simula- 
tions conducted  with  the  downwash  effects  removed  from  the  equations 
of  motion. 
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A CONTROL  THEORETIC  STUDY  OF  THE  NUCLEAR 


VULNERABILITY  OF  THE  AIR  REFUELING  TASK 
I.  Introduction 

There  is  much  current  Interest  in  determining  the  ability  of  the 
pilot  of  a B-52  aircraft  performing  a strategic  bombing  mission  to 
suffer  the  radiation  effects  of  a nuclear  encounter  and  still  continue 
on  to  mission  completion.  A typical  mission  consists  of  a takeoff, 
cllfflbout,  cruise,  air  refueling,  low  level  phase,  and  a penetration 
to  the  target.  The  air  refueling  portion  of  the  mission  is  one  of 
the  most  vulnerable  to  the  effects  of  pilot  irradiation  because  it 
requires  flying  in  close  formation  with  a tanker  aircraft  in  order 
to  receive  enough  fuel  to  complete  the  mission.  This  study  treats 
the  air  refueling  portion  of  the  B-52  mission. 

The  longitudinal  small  perturbation  equations  of  motion  of  the 
B-52H  aircraft  were  programmed  on  an  analog  computer  and  a two- 
dimensional  simulation  of  the  air  refueling  was  made.  These  linear- 
ized equations  Included  the  force  and  moment  Induced  on  the  B-52  by 
the  KC-135  tanker  downwash  flow  field.  The  pitch  stability  augmenta- 
tion system  (SAS)  and  the  engine  dynamics  also  were  modeled  on  the 
analog  computer. 

Pilot  models  were  developed  which  enabled  a closed-loop  analysis 
of  the  elevator  and  throttle  control  loops  of  the  simulated  B-52. 

The  parameters  of  these  pilot  models  were  chosen  on  the  basis  of 
human  operator  theory  and  stability  considerations.  To  this  end  It 
was  necessary  to  perform  frequency  domain  analyses  on  the  transfer 
functions  of  the  controlled  elements  and  the  pilot-aircraft  open-loop 


transfer  functions.  Although  the  analog  simulation  was  concerned 
primarily  with  Che  pilot-aircraft  systems  having  downwash  effects 
Included,  the  frequency  domain  analyses  were  performed  with  equal 
emphasis  on  the  systems  with  and  without  downwash. 

The  variables  considered  In  this  study  are  listed  below: 

1.  Pilot  lead  time, 

2.  RMS  vertical  gust  intensity, 

3.  Thrust  (throttle)  authority  limits 

4.  Radiation  dose  level 

The  pilot  lead  was  the  only  variable  In  the  theoretical  analyses, 
whereas  all  four  variables  were  considered  In  the  analog  simulation. 

The  primary  objective  of  this  study  was  to  predict  the  length  of 
time  required  to  accomplish  an  air  refueling  under  the  adverse  condi- 
tions of  air  turbulence  gusts  and  pilot  performance  degradation  due  to 
irradiation.  Determination  of  the  values  of  Tj^  which  produced  the 
lowest  values  of  predicted  refueling  time  was  a secondary  objective. 
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The  Coordinate  Systems 

Three  coordinate  systems  were  used  in  this  study.  These  are 
depicted  In  Figure  1.  The  X,Z  system  has  its  origin  at  the  boom  pod 
of  the  tanker  aircraft.  These  axes  are  aligned  with  the  local  hor- 
izontal and  vertical.  The  perturbed  displacements  of  the  B-52  center 
of  mass  along  these  axes  are  x and  z,  respectively.  Since  the  dis- 
placements of  the  B-52  boom  receptacle  are  of  primary  Interest  to  this 
problem,  these  displacements,  x'  and  z',  are  defined  as 


i ! 


X + 54  cos0 
z - 54  8in0 


(1) 

(2) 


where  54  Is  the  distance  in  feet  from  the  B-52  center  of  mass  to  the 
refueling  boom  receptacle  and  0 Is  the  perturbed  pitch  attitude.  The 
B-52  fuselage  Is  assumed  here  to  have  zero  vertical  dimension.  If  It 
is  further  assumed  that  all  perturbations  In  0 are  small,  the  above 
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The  R,n  coordinate  system  Is  derived  from  the  X,Z  system  by  the 


following  relationships: 


« /-2  ^ '2 
R » /x  +2 


n ■ tan 


-1  z' 


where  R Is  the  magnitude  of  the  radius  vector  from  the  tanker  boom  pod 
to  the  B-52  boom  receptacle  and  n is  the  angle  of  the  radius  vector 
from  the  local  horizontal.  When  the  refueling  boom  is  connected  to 
the  B-52,  R and  n represent  the  boom  length  and  boom  angle,  respectively. 

The  refueling  envelope  is  the  region  in  space  defined  by 


33.67’  -^R-  45.92' 


20°  i ^ i 40*^ 


(Ref  1:1-6) 


The  X,Z  system  of  coordinates  has  its  origin  at  the  center  of  the 
refueling  envelope.  These  axes  also  are  aligned  with  the  local  hor- 
izontal and  vertical.  Since  small  perturbations  In  6 are  assumed,  the 
B-52  stability  axes  always  remain  approximately  aligned  with  the  X,Z 
axes.  R and  2 are  the  perturbed  displacements  of  the  B-52  center  of 
mass  from  the  center  of  the  refueling  envelope.  R'  and  i'  are  the 
perturbed  displacements  of  the  boom  receptacle  from  the  center  of  the 
envelope,  where 


z'  - 2 - 540 


which  are  derived  In  a manner  Identical  to  the  derivation  of  Eqs  (5) 
and  (6). 
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The  X,Z  system  Is  related  to  the  X,Z  system  by  the  relationships 
listed  below. 


X = X + 34.47  ft 


z = z - 19.90  ft 


(13) 

(14) 


The  precontact  position,  where  the  B-52  is  stabilized  before  at- 
tempting hookup,  is  depicted  on  Figure  1 as 


X - -50  ft 


z - +10  ft 


(Ref  1:4-2) 


(15) 

(16) 


Table  I,  below,  presents  the  coordinates  of  the  refueling  envelope 
center  and  the  precontact  position  in  all  three  coordinate  systems. 

TABLE  I 

Position  Coordinates 

Coordinate  Systems 


Position 

X.Z 

R.n 

X.Z 

ft, ft 

ft, deg 

ft, ft 

Center  of  Refueling 
Envelope 

-34.47,19.90 

39.79,30 

0,0 

Precontact 

-84.47,29.90 

89.61,19.5 

-50,10 

The  Air  Refueling  Task 

The  air  refueling  phase  of  the  hypothetical  strategic  bombing  mis- 
sion was  assumed  to  take  place  approximately  five  and  one  half  hours 
after  takeoff.  The  B-52  takes  on  120,000  lbs  of  fuel  during  an  "on- 
boom"  time  of  twenty  minutes.  The  aircraft  parameters  and  flight  con- 
ditions assumed  for  that  point  in  the  mission  are  given  in  Table  II. 
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Table  II 

B>-52H  Longitudinal  Aerodynamic  Characteristics  for  the  Refueling  Task 

(Stability  Axis  System) 


0719  deg"-^  Z - -3.798  ft/sec 


Also  listed  in  that  table  are  the  B-52  stability  derivatives  which 
will  be  used  in  the  equations  of  motion  in  later  sections. 

The  air  refueling  maneuver  requires  the  pilot  to  perform  three 
types  of  procedures.  They  are: 

1.  hookup 

2 . boom  tracking 

3.  inadvertant  disconnect,  if  required 

The  air  refueling  task  begins  with  the  B-52  stabilized  in  the  pre- 
contact position,  depicted  in  Figure  1.  The  initial  hookup  procedure 
is  begun  when  the  pilot  advances  the  throttles  enough  to  establish  a 
slow  closure  rate  with  respect  to  the  tanker.  The  pilot  climbs  until 
the  B-52  boom  receptacle  is  slightly  below  the  refueling  boom,  which 
has  been  extended  into  the  refueling  envelope.  As  the  B-52  boom  re- 
ceptacle approaches  the  center  of  the  refueling  envelope  the  pilot 
reduces  power  in  order  to  arrive  at  the  center  of  the  envelope  with 
zero  closure  rate.  The  boom  operator  aboard  the  tanker  then  Inserts 
the  boom  nozzle  into  the  B-52  boom  receptacle  (Ref  1:4.2). 

The  boom  tracking  procedure  requires  the  B-52  pilot  to  keep  his 
aircraft  within  the  limits  of  the  refueling  envelope  until  the  refuel- 
ing Is  completed  or  until  the  occurrence  of  an  Inadvertant  disconnect. 

An  Inadvertant  disconnect  occurs  when  the  boom  receptacle  moves 
outside  the  limits  of  boom  travel  as  defined  by  the  refueling  envelope. 
If  refueling  is  not  completed  the  B-52  pilot  returns  his  aircraft  to  the 
precontact  position  and  stabilizes  before  attempting  another  hookup. 
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Assumptions.  Approximations,  and  Simplifications 

The  complexity  of  the  air  refueling  problem  required  that  the  task 
be  reduced  to  a simpler  form  for  this  study.  The  statements  listed  I 

below  define  the  manner  in  which  this  was  done. 

1.  The  B-52  was  treated  as  a rigid  body.  1 

2.  The  weight  of  the  B-52  remained  constant  at  450,000  lbs  , 

during  refueling.  It  was  assumed  that  maximum  gross  weight  represented 

the  worst  case  from  the  standpoint  of  pilot  control. 

3.  The  KC-135  tanker  was  treated  as  a stable  platform,  unaf- 
fected by  the  vertical  air  turbulence  gusts. 

4.  The  B-52  was  assumed  to  have  accomplished  a successful 
hookup  when  it  had  remained  for  fifteen  seconds  within  an  envelope 
defined  by 

37.79'<R<  41.79’  (17) 

25°  < n < 35°  (18) 

This  envelope  is  shown  by  dashed  lines  within  the  larger  refueling 
envelope  in  Figure  1. 

5.  In  the  Inadvertant  disconnect  procedure,  the  B-52  was 

V 

assumed  to  have  stabilized  in  the  precontact  position  when  It  had 
remained  for  fifteen  seconds  within  an  envelope  defined  by 


85 ' <.  R £ 

95' 

(19) 

15°  < n < 

25° 

(20) 

This  envelope  also  Is  defined  by  dashed  lines  In  Figure  1. 
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6.  During  the  hookup  procedure  the  B-52  closure  rate  was 
required  to  remain  between  one  and  two  feet  per  second  relative  to  the 
tanker. 

7.  The  B-52  was  assumed  to  be  affected  only  by  vertical 
air  turbulence  gusts.  Head-on  and  side  gusts  were  not  considered. 

8.  The  derivative  of  the  vertical  gust  velocity,  Wg,  was 

neglected. 

9.  Small  perturbations  in  u,w,q,  and  6 were  assumed. 

10.  The  B-52  stability  axes  were  assumed  to  be  always  ap- 
proximately aligned  with  the  local  horizontal  and  vertical. 

11.  The  engine  dynamics  were  approximated  by  a first  order 
lag  transfer  function. 

12.  The  B-52  pitch  SAS  was  simplified  to  that  shown  in 
Figure  2 by  neglecting  several  high  frequency  terms  in  the  transfer 
functions  of  its  components  (Ref  3:17). 

13.  A linear  form  of  pilot  model,  Yp  * Kp(TLjuH-l)e”^“‘^e\  was 
assumed  both  for  the  elevator  and  the  throttle  control.  This  model 
will  be  the  subject  of  discussion  in  later  sections. 

14.  The  terms  in  the  pilot  models  were  represented  in 

the  analog  simulations  by  second  order  Fade' approximations . 

15.  The  pilot  models  applied  control  inputs  directly  into 
the  pitch  SAS  and  the  engine  throttles.  Connecting  linkages  and  gear 
reductions  normally  present  between  the  cockpit  and  these  control 
units  were  disregarded.  The  gains  of  the  pilot  models  were  calculated 
on  the  basis  of  these  approximations. 

16.  The  B-52  was  assumed  to  be  trimmed  at  all  times  for 
flight  at  the  center  of  the  refueling  envelope.  As  shown  later  this 


assumption  simplified  the  dovmwash  equations. 

17.  The  B-52  was  assumed  to  be  aligned  at  all  times  on  the 
centerline  of  the  KC-135  tanker  with  all  lateral  variables  equal  to 


18.  The  hypothetical  nuclear  encounter  was  assuiaed  to  occur 
at  the  time  of  takeoff. 
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III.  Analytical  Development 

The  Pilot-Aircraft  Control  System 

A block  diagram  of  the  pilot-aircraft  control  system  is  shown  in 

Figure  2.  ¥„  is  the  transfer  function  of  the  pilot  model  for  elevator 

control  and  Yn  is  the  pilot  model  for  throttle  control. 

Ft 

The  stability  characteristics  of  the  closed-loop  control  system  can 

be  analyzed  if  Yj,  Y.  and  Y„  Y^ 

P'j  the  open-loop  pilot-aircraft 

transfer  functions,  are  derived.  Y-  and  Y.  , the  transfer  functions 

E 

of  the  elevator  and  throttle  controlled  elements,  take  the  following 
forms  for  the  air  refueling  task: 


(21) 

(22) 


where  5^  represents  the  pilot's  input  to  the  pitch  SAS  and  S is  the 

*^E  Pt 

pilot's  input  to  the  engine  throttle. 

The  system  shown  in  Figure  2 can  be  modeled  in  terms  of  pilot 
models  operating  on  controlled  elements  as  shown  by  Figure  3.  The 
succeeding  sections  will  lead  up  to  the  derivation  of  the  terms  of 
the  two  controlled  elements. 


Engine  Dynamics 

At  the  refueling  altitude  of  32,700  ft  the  thrust  required  to 
maintain  equilibrium  flight  is  calculated  from  the  equilibrium  drag 
coefficient,  Cq^,  as  shown  below. 
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Reduced  Block  Diagram  of  Pilot-Aircraft  Control  System. 


Teq  - Drag  » = (.0278)  (198)  (4000) 


22017.6  lbs 


Equilibrium  thrust  per  engine  = 22017.6/8  = 2752.2  lbs 


The  above  value  of  equilibrium  thrust  defines  an  operating  point  about 
which  positive  or  negative  changes,  6,p,  are  commanded  by  the  pilot. 

The  thrust  time  response  of  the  TF33-P-3  engine  of  the  B-52H  for 
a step  throttle  movement  from  idle  to  military  rated  thrust  (MRT)  is 
plotted  in  Reference  4.  That  portion  of  the  curve  representing  the 
time  response  of  to  a step  throttle  movement,  , from  the  equilib- 

rium thrust  setting  to  MRT  is  reproduced  in  Figure  4.  By  trial  and 
error  it  was  determined  that  the  time  response  for  this  particular 
case  can  be  approximated  by 


firCt)  “ « (l-e~^‘=)  = 926 d-e'^*^) 

T p^ 


which  is  also  plotted  on  Figure  4 for  comparison.  Taking  Laplace 
transforms  of  both  sides  of  the  above  expression  yields  the  transfer 
function 


.25  s + 1 


This  transfer  function  is  used  in  the  analog  simulations  to  represent 
the  engine  dynamics.  It  is  by  no  means  meant  to  provide  a close  approx- 
imation to  the  engine  response  in  the  general  case.  On  the  contrary, 
the  approximation  was  based  on  the  response  to  only  one  of  the  many 
possible  throttle  inputs — that  of  a step  input  from  equilibrium  to 
MRT  thrust.  A step  throttle  input  from  idle  thrust  to  the  equilibrium 
thrust  setting,  for  example,  will  result  in  a much  slower  engine 
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Figure  4.‘  Comparison  of  B-52H  Thrust  Time  Response  with  the  Time 
Response  of  the  Approximation  (Ref  4 >142). 


I 


response  Chan  given  by  the  approximation.  However,  It  Is  assumed 
that  Che  approximation  of  Che  engine  response  In  the  neighborhood  of 
the  equilibrium  thrust  setting  will  also  be  a valid  approximation  to 
the  engine  response  In  a statistical  sense;  l.e.,  the  average  response 
experienced  over  many  simulations  of  the  air  refueling. 


The  Downwash  Model 

The  effects  of  the  KC-135  flow  field  on  the  B-52  were  derived 
from  graphical  data  in  the  form  of  induced  changes  In  the  B-52  equilib- 
rium lift  and  moment  coefficients,  and  (Ref  5:2-79,4-95). 

This  data  is  given  as  a function  of  the  variable  z.  A linear  fit  is 
made  to  this  data  as  shown  in  Appendix  A.  The  resultant  equations  for 
the  induced  changes,  converted  to  the  variable  z,  are 

ACr,  - -.077  + .OOlUz  (26) 


- .068  - .00102Z  ( 

This  study  assumed  that  the  effects  of  tanker  downwash  are 
trimmed  out  when  the  B-52  Is  flying  at  the  center  of  the  refueling 
envelope.  In  that  case  the  conditions 


z-o  " 0 


“1  £-o 


must  hold  since  the  aircraft  must  have  the  same  total  lift  ( > weight) 
and  moment  (’>  0)  as  In  the  case  with  no  downwash.  Therefore,  when 
the  aircraft  Is  trimmed  for  flight  at  the  center  of  the  refueling 
envelope  the  constant  terms  drop  out  of  the  downwash  equations, 

A 

reducing  them  to  simple  linear  terms  In  z: 
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- .00114z  (30) 

“ -.00102Z  (31) 

The  changes  In  lift  and  moment  caused  by  the  tanker  downwash  are 

AL  » ACl^  S (32) 

AM  = AC^^  qi  S c (33) 


where  is  the  dynamic  pressure,  S is  the  wing  area,  and  "c  is  the 
mean  geometric  chord  of  the  wing. 

In  order  to  trim  out  the  induced  effects  of  the  tanker  flow 
field,  the  B-52  pilot  must  change  the  steady  state  angle  of  attack 
by  changing  the  horizontal  stabilizer  angle  by  an  Increment 
dstab*  ^ stability  axis  system  Is  thus  defined.  The  steady  state 
angle  of  attack  which  Is  required  for  flight  at  the  center  of  the 
refueling  envelope  Is  .68  deg,  as  shown  by  the  calculations  In  Ap- 
pendix B.  This  value  of  was  assumed  to  be  small  enough  to  be 
treated  as  equal  to  zero  In  the  ensuing  derivation  of  the  B-52H 
equations  of  motion. 

B-52H  Equations  of  Motion 

The  B-52H  equations  of  motion,  using  the  dimensional  form  of  the 
stability  derivatives,  are  presented  below.  These  equations  Include 
the  vertical  gust  disturbance,  Wg,  and  the  lift  and  moment  terms 
Induced  by  the  KC-135  downwash. 

u - X„  u + X„(w-Wg)-g6  + 6^  (34) 
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(35) 


Zo.  qjS 

4-MuU+M*,(w-w„)-«lqq+^  w +MA(w-w_)+M  6r+M.  (36) 

Ui  * g 6E  *T  lyy  ■" 

(Ref  6:9.83) 

Neglecting  the  Wg  terms,  rearranging  the  equations,  and  substituting 


In  values  for  the  coefficients  as  listed  In  Table  II,  yields 

u - -.00A465U  +.02548W  - 32.26  - .025A8wg  + .00007156^  (37) 

* - -.09143U  - .4117w  + 701. 37q  + .4ll7wg  - 8.1436g  + (38) 

q - .000013U  - .002676W  - 1.00505q  +.002676Wg  - .626555g 

+ 3.411xl0"^6j  + Bz  (39) 

0 - q (40) 


where  the  coefficients  of  the  downwash  variable,  z,  are: 

-qis 

A - —4-  (.00114)  = -.06471  (41) 

00102)  = -.002269  (42) 

“ ~T 

lyy 

which  are  derived  from  Table  II  and  Eqs  (30)  and  (31) . 

By  introducing  the  equation, 

I - w - Ui  0 (43) 

and  by  letting  u,  w,  q,  8,  z be  the  state  variables  and  dg,  be 
the  control  variables  of  the  system,  the  equations  of  inotlon  can  be 
put  Into  matrix  format  as  shown  below. 
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^1 ' 

-.004465 

.02548 

0 

-32.2 

0 

u 

-.09143 

-.4117 

701.37 

0 

A 

w 

q 

► = 

.000013 

-.002676 

-1.00505 

0 

B 

< 

q 

e 

0 

0 

1 

0 

0 

0 

, 

0 

1 

0 

-706 

0. 

■-.02548  * 

0 .0000715  ■ 

.4117 

-8.143  0 

.002676 

Wg  + 

-.62655  3.411x10"^ 

0 

0 0 

0 

0 0 

(44: 


The  characteristic  equation  of  this  system  Is 

C.E.  - s^  + 1.4212s^  + (2.2994-A)s^  + (. 012756-1. 00953A+4.6B)s^ 

+ (.0080505-1. 89374A+290.6807B)s  + (-.00862A-.001525B)  (45 

The  eigenvalues  and  eigenvectors  of  the  system  without  the  down- 
wash  effects  were  calculated  by  letting  A and  B go  zero,  reducing  the 
system  to  fourth  order.  The  resultant  characteristic  equation  Is 

C.E.jjoDW  * s^l.4212s3+2.2994s2+.012756s+. 0080505  (46: 

Substituting  the  numerical  values  for  A and  B Into  Eq  (45) 
yields  the  characteristic  equation 

C.E.jjj^  - 8^+1. 4212s V2.3641s^+.067645s^-.52896s+. 0005612 


which  represents  the  system  with  downwash.  The  eigenvalues  and  eigen' 
vectors  for  this  system  also  were  calculated.  A comparison  of  the 
elgenanalysls  for  the  B-52  with  and  without  downwash  Is  presented  In 


Table  III.  The  downwash  terms  introduce  two  poles  In  the  right  half 


8-plane,  thus  causing  the  aircraft  to  be  unstable.  Additionally,  the 
phugold  oscillatory  roots  are  eliminated. 


Derivation  of  the  Controlled  Elements 

The  controlled  elements  Ycg  and  have  been  defined  as 


_g_  = (£_  - 54  _3_ 

*^Pe  '^E  "^Pe  '*Pe 


CT  5p^  6^  6p,j  s 6p^ 


where  -x — is  the  transfer  function  of  the  engine  dynamics  and 


Y(._^  represent  the  vertical  and  horizontal  displacements  of  the  B-52 
boom  receptacle  due  to  the  respective  pilot  Inputs  6p^  and  6^^.  The 
throttle  controlled  element  can  be  derived  in  a straightforward  manner 
from  the  equations  of  motion,  but  the  derivation  of  the  elevator  con- 
trolled element  requires  a complicated  block  diagram  reduction  of  the 
pitch  SAS  feedback  control  loop.  Details  of  the  derivation  of  the  con- 
trolled elements  are  shown  in  Appendix  C. 

Table  IV  presents  a summary  of  the  controlled  elements  in  factored 
frequency  domain  form.  Inspection  of  these  controlled  elements  re- 
veals that  the  downwash  terms  have  caused  them  to  become  unstable  by 
Introducing  two  positive  real  poles.  In  addition  to  eliminating  the 
phugold  oscillation,  the  tanker  downwash  changes  the  short-period 
response  of  the  elevator  controlled  element  by  decreasing  Its  damping 
while  increasing  Its  damped  frequency. 


Table  HI 

Eigenanalysis  of  B-52H  Characteristic  Equation 


Without  Downwash 


With  Downwash 


EIGVAL(I) 

-.7089110  1.337319 


EIGVAL  (1) 

-.6379154  1.411451 


EIGVEC  (1 


.1097337D-01 

1.000000 

-.422327UD-03 

.1244091D-02 


.5121876D-02 

.0 

.1907391D-02 

.3436895D-03 


EIGVAL  (2) 


-.7089110 

EIGV 

.1097337D-01 

1.00000C 

-.422327UD-03 

.1244091D-02 


-1.337319 

EC(2) 

-.5121876D-02 

.0 

-.19073910-02 

.34368950-03 


EIGVEC  (1) 


. 10628420-01 
1.O00000 
•.  31046890-03 
.12598580-02 
.1157477 


.59193350-02 

.0 

.20011390-02 
-.34943750-03 
-.  1306298 


-.6379154 


EIGVAL (2) 

-1.411451 


EIGVEC  (2) 


. 10628420-01 
1.00000C 
-. 31046890-03 
.12598560-02 
.1157477 


-.59193350-02 

.0 

-.20011390-02 
.34943750-03 
. 13C6298 


EIGVAL (3) 

-.16964630-02  .59255150-01 


EIGVEC  (3) 


1.000000 

-.36199530-01 

.10928600-03 

-.11462430-03 


.0 

-.10443390-02 

-.36688060-05 

-.18410480-02 


EIGVAL (3) 
.10612460-02  .0 

EIGVEC  (3) 
.9193964  .0 

-.2516848  .0 

-.37881010-06  .0 

-.35694850-03  .0 

.3022661  .0 


EIGVAL  (4) 

-.16964630-02  -.59255150-01 


EIGVEC (4) 


1.000000 

-.36199530-01 

.10928600-03 

-.11462430-03 


.0 

.10443390-02 

.36688060-05 

.18410480-02 


EIGVAL (4) 
-.5484181  .0 

EIGVEC (4) 
.96302910-02  .0 

-.6662926  .0 

.19992920-03  .0 

-.36455620-03  .0 

.7456281  .0 


EIGVAL (5) 
.4019726  .0 


EIGVEC  (5) 


.63422050-01  .0 
-.4873993  .0 
-.47682490-03  .0 
-.11862120-02  .0 

.8708719  .0 


Note:  Elements  of  eigenvectors  correspond  to  elements  of  state  vector 
of  Eq  (44). 


Table  IV 


Transfer  Functions  of  the  Controlled  Elements 


NoDW 

) 

T NoDW 

E DW 

Yct) 

T DW 

Numerator 

K 

-7.766 

.02034 

-112.237 

-.06857 

Tl 

6.667 

- 

-2500 

1.82 

T2 

4.0 

- 

- 

-2.453 

T3 

-2500 

- 

- 

- 

“ni 

2.544 

- 

2.557 

- 

5l 

.058 

- 

.058 

- 

Denominator 

Ts 

1 

- 

- 

1 

^1 

4.5 

.25 

1.585 

.25 

^2 

6.0 

- 

-2.57 

1.824 

T3 

- 

- 

-952.38 

-2.488 

T4 

- 

- 

-943.4 

“1 

1.514 

.0593 

1.471 

- 

00 

• 

.028 

.62 

- 

% 

.0593 

- 

- 

- 

«2 

.028 

- 

- 

- 

General  form  of  numerator  or  denominator: 

r.  - K(T,J»)  jSj  (TiJ^l)  ^ Jurt-1) 

k K 


r 


I : 
[■ 
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The  Gust  Model 

As  previously  stated,  the  B-52  is  assumed  to  be  affected  only  by 
vertical  gusts.  The  random  vertical  gust  model  used  in  this  study  is 
defined  by  the  following  transfer  function,  which  is  based  on  the 
Dryden  scales  for  clear  air  turbulence. 


where 


'w 


2.85  (A.291aH-l) 
(2.48ju>4-l)2 


(50) 


Wg  * random  vertical  gust  velocity 


o„  * rms  value  of  random  vertical  gust  velocity 


C rms  value  of  Gaussian  random  noise  input 

(Ref  7:459) 

The  bandwidth  of  the  random  noise,  was  1.5  Hz-  Calculation  of  the 
terms  in  this  gust  model  is  shown  in  Appendix  D. 

Bode  magnitude  and  phase  angle  plots  of  the  gust  transfer  function 
are  shown  in  Figure  D-1.  It  is  seen  that  the  log  magnitude  of  the 
transfer  function  is  at  its  half-power  point,  or  3 db  down,  at  a 
frequency  of  approximately  .82  rad/sec.  This  frequency  was  considered 
to  be  the  highest  frequency  at  which  effective  gusts  could  be  generated 
by  the  model. 

The  Pilot  Model 

The  form  of  the  pilot  model  which  was  chosen  for  this  air  refueling 
simulation  is  the  linear  describing  function 


tp  - Kp  (Tl  jori-  1)  e'J“^e 


(51) 


t I 

t 1 

I j 
i { 


( ; 
5 .1 


24 


I 


where  Kp  = pilot  gain  1 

Tl  “ lead  time  constant,  sec  j 

Tg  * effective  time  delay,  sec 

This  model  has  been  used  with  success  in  previous  work  with  compensatory 
control  systems  (Ref  8:4).  The  parameters  Kp,  Tj^  and  Tg  must  be 

chosen  as  those  which  best  fit  each  particular  control  problem  under  \ 

consideration  whether  it  be  the  air  refueling  with  downwash  or  without  ! 

downwash . > 

The  pilot's  time  delay,  Tg,  was  chosen  as  0.3  seconds  for  all 
cases.  This  value  of  time  delay  has  been  used  in  many  other  pilot 
modeling  applications  and  has  been  verified  in  many  human  operator 
tracking  experiments.  It  may  be  possible  that  for  the  air  refueling 
task  the  human  pilot  could  work  at  a lower  value  of  time  delay,  but 
there  is  currently  no  data  in  existence  to  support  that  theory. 

Since  the  controlled  elements  under  consideration  here  are  un- 
stable (with  downwash)  and  are  more  complicated  than  those  investigated 
In  many  of  the  previous  human  operator  tracking  experiments,  their 
forms  did  not  suggest  appropriate  values  for  T^.  It  has  been  shown 
that  in  an  actual  flight  environment  a pilot  may  generate  a value  of 
Tx.  ss  high  as  eight  seconds  (Ref  8:4).  For  this  reason,  the  value  of 
Tj^  was  varied  over  a range  of  from  one  to  eight  seconds  to  determine 
which  values  provided  the  best  performance  in  the  air  refueling  task. 

There  are  indications  from  previous  work  in  the  literature  that  a 
highly  trained  pilot  would  be  operating  with  the  optimum  Tj.,  if  one 
exists. 


The  pilot  gains,  Kp,  were  calculated  Individually  for  each  value 
of  Stability  considerations  were  the  primary  bases  for  determining 

values  of  Kp,  as  shown  in  the  following  section. 

It  has  been  found  that  the  human  pilot  adapts  the  form  of  his 
equalizing  characteristics  to  achieve  stable  control,  good  low  fre- 
quency closed  loop  system  response  to  the  forcing  function,  and 
maintenance  of  high  frequency  system  stability  (Ref  9:18-19).  These 
considerations  were  important  in  choosing  the  parameters  of  the  pilot 
models  used  in  this  study. 


Analysis  of  the  Pilot-Aircraft  System 


In  this  section  the  open-loop  transfer  functions  TpY^,  are 
analyzed  in  the  frequency  domain  in  order  to  determine  the  stability 
characteristics  of  the  corresponding  closed-loop  systems.  The  ob- 
jective of  this  analysis  is  to  find  values  of  pilot  gain  , Kp,  which 
are  most  suitable  with  respect  to  stability  and  system  following 
characteristics  as  T^  is  varied.  The  Nyqulst  stability  criterion  is 
the  basis  for  these  stability  analyses  (Ref  10:272-287), 

The  pilot-aircraft  open-loop  transfer  functions  for  the  case 
without  downwash  are  analyzed  first.  These  transfer  functions  are 
shown  below  for  the  elevator  and  the  throttle  control  loops. 


^o„Yc„)  - Ko  (TTj<iH-l)e"^“-3  HfcK^  1 (52) 

uo  DW  ^ L ‘'E  No  Dwj 

Wcn)  - KpT(TLJ<^l>e"J“*3  fYr  ) 1 (53) 

‘=T''no  DW  L |_  no  DwJ 


Y.  ) and  Yc„)  are  given  by  Table  IV.  Bode  plots  and  polar 

Vno  DW  ‘'T^No  DW 

plots  of  these  open-loop  transfer  functions  are  provided  in  Appendix  E 


T-'"  ■ w 


for  values  of  from  1 to  8.  Applying  the  Nyqulst  criterion  to  the 
polar  plots  showed  that  the  corresponding  closed-loop  pilot-aircraft 
systems  are  stable  for  values  of  pilot  gain  between  zero  and  certain 
maximum  limits,  above  which  the  systems  are  unstable.  These  maximum 
limits  are  given  in  Table  V, 

Values  of  Kp^  and 

Rule  of  Thumb"  of  frequency  response  synthesis,  stated  below: 


K_  were  calculated  by  applying  the  "Primary 
Ft 


"Find  or  create  a fair  stretch  of 
-20db/decade  slope  for  the  amplitude 
ratio,  and  then  make  it  the  crossover 
region  by  putting  the  0 db  line 
through  it"  (Ref  9:27). 


Additional  considerations  when  using  this  method,  in  order  of  decreasing 
Importance,  are 


1.  The  system  must  have  a positive  phase  margin. 

2.  For  a "good"  feedback  control  system  the  pilot's  describing 
function,  Yp,  must  be  adjusted  so  that  the  crossover  frequency,  u>^, 
exceeds  the  highest  Important  frequency  of  the  forcing  function  (Ref 
9:27).  The  forcing  function  in  this  problem  is  the  vertical  gust,  Wg. 
It  has  already  been  shown  that  the  -3  db  frequency  of  the  gust  model 

Is  .82  rad/sec.  If  this  frequency  is  taken  to  be  the  highest  Important 
frequency  of  the  input  then  the  crossover  frequencies  of  the  open-loop 
transfer  functions  should  always  exceed  .82  rad/sec. 

After  using  the  Primary  Rule  of  Thumb  to  define  the  crossover 
region,  the  values  of  Kn  were  calculated  for  a crossover  frequency 
of  .82  rad/sec.  These  are  given  in  Table  V.  The  systems  are  shown  to 
be  unstable  for  T^'s  of  1.5  and  1 sec.  Those  two  cases  would  require 
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Table  V 


Summary  of  Pilot-Aircraft  System  Analysis 
^PE  ^ No  DW  ^Pt 


Y Y 1 
Pt  VNo  DW 


1 

Tl 

<^c 

s 

“c 

4>m 

8 

.001698 

.82  r/s 

23  deg 

753 

.45  r/s 

61  deg 

(0-. 002927) 

(0-6469) 

i 

7 

.001905 

.82 

22 

933 

.49 

59 

} 

1 

(0-. 003302) 

(0-7338) 

i 

i 6 

.002238 

.82 

20 

1124 

.51 

56 

« ; 

(0-. 003706) 

(0-8477) 

I 

5 

.002660 

.82 

18 

1520 

.58 

53 

F 

(0-. 004222) 

(0-10034) 

f 

4 

.003283 

.82 

15 

2031 

.62 

49 

1 

(0-. 004953) 

(0-12291) 

3 

.004240 

.82 

9 

2991 

.70 

43 

I 

. 

(0-. 005723) 

(0-15859) 

\ 

i 

t 2 

.005752 

.82 

1 

5029 

.83 

33 

1 

! 

(0-. 005934) 

(0-22347) 

i 

1.5 

.007076 

.82 

-7 

6848 

.90 

26 

(0-. 004250) 

(0-26720) 

b 

1 

.008756 

.82 

-20 

9111 

.93 

14 

(0-. 000303) 

(0-30921) 

■•J 

’i 

Note: 

Values  In  parentheses 

are  the  minimum  and 

maximum  values 

of 

stable  K . 
P 


(Table  V continued  on  next  page) 
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Table  V 
(contd) 


^Pr  ^DW  ^PT  ^DW 


Tl  KpE 

“c 

KpT 

liiC 

8 .003531 

(.003406-. 003667) 

1.24  r/s 

2 deg  794 
(1.9  - 

CD  ) 

.47  r/s 

61  deg 

7 .003731 

(.003406-. 004125) 

1.19 

6 

939 
(2.2  - 

CD  ) 

.49 

59 

6 .003953 

(.003404-. 004714) 

1.12 

10 

1183 
(2.7  - 

CD  ) 

.53 

56 

5 .004195 

(.003403-. 005466) 

1.03 

16 

1526 
(3.4  - 

CO  ) 

.58 

53 

4 .004467 

(.003402-. 006503) 

.90 

23 

2037 
(4.4  - 

CD  ) 

.62 

49 

3 .004695 

(.003402-. 007816) 

.70 

30 

2996 
(5.9  - 

CD  ) 

.70 

43 

2 .004924 

(.003402-. 008910) 

.47 

26 

5033 
(9.8  - 

CD  ) 

.83 

33 

1.5  .004925 

(.003403-. 008910) 

.39 

18 

6848 
(14.9  - 

CD  ) 

.90 

26 

1 .004929 

(.003407-. 008910) 

.34 

10 

9115 
(30.4  - 

CD  ) 

.93 

14 

NbCe:  Values  In  parentheses  are  the  mlnlnnim  and  maxlmuni  values  of 

stable  Kp. 


Wq  to  be  decreased  to  a value  less  than  .82  rad/sec  In  order  to  achieve 
stability.  The  values  of  Kp^  were  calculated  at  the  to^’s  which  pro- 
vided the  greatest  phase  margins  in  the  crossover  regions.  These 
values  also  are  presented  In  Table  V. 

The  pilot-aircraft  open-loop  transfer  functions  Including  the 
tanker  downwash  are 

''pe  Wd«  ■ ffPE  >Dw]  »4) 

*Pl  |Voj  (55) 

where  Ycg  and  can  be  derived  from  Table  IV.  Bode  and  polar 

plots  of  these  transfer  functions  are  presented  in  Figures  F-1  through 
F-18. 


Nyquist  analysis  of  the  polar  plots  of  YPE  ) showed  that 
the  systems  with  downwash  are  conditionally  stable.  There  exists  a 
region  of  stability  within  each  polar  plot  corresponding  to  a particu- 
lar range  of  Kp^.  Values  of  KpE  for  which  the  -1+jO  point  lies  outside 
this  specific  region  cause  the  pilot-aircraft  system  to  be  unstable. 
These  regions  of  stability  are  marked  with  cross-hatching  on  the  polar 
plots.  It  can  be  seen  that  the  stable  region  is  quite  small  for 
Tl  ■ 8 and  increases  in  size  as  the  value  of  T^  decreases.  The  values 
of  Kpg  which  placed  the  -1+jO  points  at  the  graphical  centers  of  the 
stable  regions  were  calculated.  These  values,  given  in  Table  V,  were 
used  in  the  analog  simulations.  Table  V also  lists  the  range  of  stable 
Kpg  for  each  value  of  Tj^. 


Applying  the  Nyquist  stability  criterion  to  the  polar  plots  of 


Yp^  Yc^  ) Qy  in  Appendix  F showed  that  these  transfer  functions  do  not 
exhibit  the  conditionally  stable  characteristic.  They  are  unstable 
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only  If  the  pilot  gain  falls  below  some  minimum  value.  Lacking 


( 

» 


I ' 
1 


I 

1 

f 

[ 


stability  considerations  other  than  to  keep  Kp^  above  this  minimum 
value,  the  Primary  Rule  of  Thumb  was  employed  to  calculate  the  values 
of  Kp^.  As  In  the  case  of  Yp^  Yj.^  ) the  values  chosen  for 

KpT  were  those  which  gave  the  maximum  phase  margins  In  the  crossover 
regions.  Table  V presents  these  gains  In  addition  to  the  minimum 
values  of  Kp^  which  are  stable. 

The  pilot  gains  calculated  In  this  section  apply  only  to  the  model 
of  the  non-lrradlated  pilot.  The  following  section  Introduces  factors 
which  can  be  applied  to  these  gains  to  simulate  performance  decrements 
due  to  Irradiation. 

Effects  of  Irradiation  on  the  Human  Pilot 

It  has  been  shown  that  a pilot  who  has  suffered  a lethal  dose  of 
radiation  Is  likely  to  suffer  a decrement  In  performance  which  varies 
as  a function  of  time  after  the  dose  Is  received.  Data  which  support 
this  conclusion  were  obtained  In  laboratory  experiments  with  trained 
monkeys  performing  flight-simulating  tasks.  The  clinical  response  to 
radiation  exhibited  by  monkeys  is  essentially  the  same  as  that  of 
humans. 

Generally  the  data  showed  that  the  monkeys  suffered  substantial 
degradations  In  performance  shortly  after  irradiation,  followed  by 
periods  of  partial  recovery  lasting  for  several  hours  (Ref  11:946). 

It  is  assumed  here  that  a decrement  In  pilot  performance  due  to 
a nuclear  encounter  can  be  simulated  by  reductions  In  the  gains 
Kpg  and  Kp^.  It  is  further  assumed  that  such  a nuclear  encounter 
occurs  at  the  time  of  takeoff.  Figure  5 illustrates  the  variation  of 
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pilot  gain  factor  as  a function  of  time  after  exposure  for  four  levels 
of  radiation  dose  (Ref  12).  The  pilot  gains  of  Table  V are  multiplied 
by  the  appropriate  gain  factor  from  this  figure  to  simulate  pilot 
Irradiation.  Table  VI  lists  the  pilot  gain  factors  which  are  applica- 
ble to  the  beginning  of  the  refueling  simulation  occurring  at  five  and 
one-half  hours  after  radiation  exposure. 


Figure  5.  Pilot  Gain  Factor  versus  Time  (Ref  12). 


Table  VI 


1 Pilot  Gain  Decrements 

Due  to  Irradiation 

Radiation  Dose 

Pilot  Gain  Factor 

(rads) 

(t  ■ si  hours)  1 

j 250 

0.96  1 

500 

0.93 

' 750 

0.87 

1000 

0.83  . 

(Ref  12)  1 
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IV.  Analog  Simulation 

Apparatus 

The  analog  simulation  was  patched  on  an  Electronics  Associates, 
Inc.  (EAI)  Model  680  analog  computer.  This  computer  was  Interfaced 
with  a Digital  Equipment  Corporation  (DEC)  PDP-12  digital  computer  by 
means  of  an  analog-to -digital  converter.  The  PDP-12  was  equipped  with 
an  RK05  disk  drive  and  an  LP08  high  speed  line  printer.  Strip  chart 
data  was  collected  by  means  of  an  EAI  8875  recorder. 

A Hewlett-Packard  HOI-3722A  Noise  Generator  supplied  the  random 
Gaussian  noise  to  the  gust  filter. 


General 

Except  where  otherwise  noted  all  experimental  work  presented  In 
this  chapter  Involved  the  simulated  B-52  as  affected  by  the  downwash 
field  of  the  KC-135  tanker. 

The  pilot-aircraft  control  system  was  patched  on  the  EAI  680  analog 
computer  according  to  the  diagram  shown  In  Figure  G-1.  Also  shown  on 
this  patching  diagram  Is  the  error  detection  and  logic  circuitry. 
Potentiometer  settings  are  listed  In  Table  X.  The  pure  time  delay  term 
In -the  pilot  model  was  represented  by  the  second  order  Fade''  approxima- 
tion. 


2 2 

-T8  . 1 - .5ts  + .083x  s 
* 1 + .5ts  + .083x^3^ 


(56) 


As  shown  earlier  in  this  report  the  frequency  domain  analyses  were 
performed  on  the  pilot-aircraft  open-loop  transfer  functions  over  a 
range  of  Tj^'s  from  1 to  8.  For  the  analog  simulation  a narrowing  of 
this  range  by  elimination  of  undesirable  Tj,'s  seemed  feasible. 


33 


The  critical  and  demanding  task  of  accomplishing  a hookup  pro- 
cedure was  chosen  as  the  basis  on  which  to  evaluate  the  performance  of 
the  pilot  model  at  the  various  values  of  t^. 

Observation  of  Hookup  Procedure 

The  characteristics  of  the  pilot-aircraft  system  accomplishing 
a hookup  procedure  in  smooth  air  (wg  = 0)  were  evaluated  by  considering 
the  following  factors: 

1.  Maximum  closure  rate  experienced  during  hookup 

2.  Hookup  flight  path 

3.  Amount  of  negative  thrust  required  to  drive  the  closure  rate 
to  zero  at  the  center  of  the  refueling  envelope 

4.  Time  required  to  accomplish  the  hookup 

To  preclude  the  danger  of  overrunning  the  tanker  during  an  at- 
tempted hookup  maneuver,  the  B-52  must  be  limited  to  a very  slow 
closure  rate.  This  restriction  was  accounted  for  in  the  analog  sim- 
ulation by  limiting  the  forward  (positive)  thrust  authority  of  the 
pilot  model  to  500  lbs.  With  this  limit  imposed  on  the  system  the 
no-gust  closure  rates  always  remained  less  than  2 ft/sec  for  all 
values  of  Tj.. 

A "good"  hookup  flight  path  was  defined  as  one  in  which  the  B-52 
gained  most  of  its  required  altitude,  with  minimum  overshoot,  before 
entering  the  refueling  envelope.  It  was  assumed  that  such  an  approach 
would  offer  the  B-52  pilot  better  visual  references  on  the  tanker 
than  would  be  possible  In  a "low"  approach. 

Negative  thrust  authority  was  limited  to  10,000  lbs  during  the 
hookup  runs.  In  other  words,  the  pilot  model  was  able  to  reduce 
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thrust  as  much  as  10,000  lbs  below  the  equilibrium  thrust  setting  In 
order  to  stop  the  closure  rate.  As  would  be  expected,  as  de^- 
creased  the  model  required  Increased  amounts  of  negative  thrust  to 
stop  the  closure  rate  due  to  Its  lesser  ability  to  "anticipate"  ar- 
rival at  the  contact  position.  Negative  thrust  requirements  more  than 
5000  lb  were  considered  excessive. 

The  time  required  for  hookup  was  that  time  required  to  enter  and 
remain  for  15  seconds  within  the  small  hookup  envelope  depleted  by 
dashed  lines  In  Figure  1. 

Figure  6 presents  the  hookup  flight  path  for  each  Tl*  These  are 
the  trajectories  of  the  B-52  boom  receptacle  as  plotted  from  the  analog 
computer.  Table  VII  gives  a comparison  of  the  maximum  closure  rates, 
negative  thrust  requirements,  and  time  to  hookup.  Based  on  the  above 
stated  criteria,  T]^'s  equal  to  3,  4,  5,  and  6 seemed  to  offer  the  best 
compromises  between  the  factors.  These  four  values  of  Tjj  were  selected 
for  further  experimentation. 

The  Simulation  Variables 

Table  VIII  lists  the  four  variables  and  the  ranges  over  which 
they  were  varied  in  the  simulations. 

The  rms  gust  Intensity,  a„,  was  varied  from  0 to  6 ft/sec  as 
shown.  The  actual  random  gust  velocities,  Wg,  produced  by  the  gust 
filter  achieved  peaks  of  roughly  three  times  the  rms  values.  Thus, 
for  o„  ■ 6 ft/sec,  the  actual  gusts  affecting  the  B-52  included  peak 
velocities  in  the  neighborhood  of  18  ft/sec. 

It  became  evident  in  the  initial  analog  runs  that  the  time 
required  to  accomplish  the  simulated  air  refueling  was  strongly 
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Table  VII 


Table  VIII 


Simulation  Variables  and  Their  Values 


Variable 

Units 

Values  Used  in  Analog  Simulations 

Tl 

sec 

3,  A,  5,  6 

ft/sec  (rms) 

0,  1,  2,  3,  A,  5,  6 

Radiation 

Dose 

rads 

0,  500,  1000  (simulated  via  pilot  gain 
decrements) 

Thrust 

lbs 

Thrust  Authority  A*,  Thrust  Authority  B** 

*Thrust  Authority  A;  1.  Hookup:  -5000  to  500  (2500) 

2.  Boom  Track:  -5000  to  2500 

3.  Inadvertant  Disconnect:  -5000  to  2500 


**Thrust  Authority  B:  1.  Hookup:  -10000  to  500  (5000) 

2.  Boom  Track;  -10000  to  5000 

3.  Inadvertant  Disconnect:  -10000  to  5000 

Values  In  parentheses  are  the  limits  of  the  graduated  forward  thrust 
for  hookup  as  explained  in  the  text. 
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Influenced  by  the  amount  of  throttle  authority  allowed  the  pilot 

model.  The  thrust  is,  of  course,  limited  by  the  capabilities  of  the 

engines  to  a range  of  approximately  -20000  lbs  to  +13000  lbs.  The  ! 

negative  limit  is  based  on  the  fact  that  not  all  of  the  22000  lbs  of 

equilibrium  thrust  can  be  eliminated  by  pulling  back  on  the  throttles; 

a small  amount  (approximately  2000  lbs)  remains  after  all  eight  ^ 

I 

throttles  are  pulled  back  to  the  idle  position.  The  positive  limit  | 

reflects  the  additional  thrust  which  is  available  as  all  eight  engine  | 

throttles  are  advanced  from  the  equilibrium  setting  (1.48  EPR)  to  j 

MRT  (2.0  EPR)  (Ref  4:104,  116).  ] 

The  decision  was  made  to  use  more  restrictive  thrust  limits  in 

the  analog  simulations.  This  decision  was  based  on  the  following  | 

i. 

thoughts : i 

I 

1.  The  B-52  may  be  flying  the  mission  with  one  or  more  engines 
Inoperative . 

2.  The  pilot  may  not  be  working  all  eight  throttles  during  the 
refueling  operation. 

3.  Adverse  engine  behavior  such  as  compressor  stalling  often 
results  from  varying  the  engine  throttles  rapidly  over  their  full 
ranges. 

Two  different  ranges  of  thrust  authority  were  treated  In  the 
analog  simulations.  The  boom  track  and  Inadvertant  disconnect  pro- 
cedures were  simulated  using  thrust  ranges  of  -5000  to  2500  lbs 
(Thrust  Authority  A)  and  -10000  to  5000  lbs  (Thrust  Authority  B).  It 
will  be  noted  that  the  negative  thrust  limits  were  two  times  as  large 
as  the  positive  limits.  Aside  from  the  fact  that  there  simply  Is  more 
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thrust  available  on  the  negative  side  of  equilibrium  than  on  the 
positive  side,  these  limits  were  established  mainly  to  fulfill  the 
desire  to  make  the  simulation  more  realistic  by  i^nposing  a passive 
type  of  safety  factor.  Since  forward  movement  of  the  B-52  is  the 
more  dangerous,  due  to  the  collision  hazard,  it  was  assumed  that  the 
pilot  would  need  more  thrust  authority  in  the  negative  direction  in 
order  to  more  quickly  stop  forward  closure  rates. 

While  the  thrust  limitations  were  easily  applied  to  the  boom 
track  and  inadvertant  disconnect  procedures,  the  simulation  of  the 
hookup  procedure  presented  some  unique  problems.  As  stated  earlier, 
the  forward  (positive)  thrust  for  hookup  was  limited  to  500  lbs  in 
order  to  keep  the  no-gust  closure  rate  at  a value  less  than  2 ft/sec. 

At  the  same  time,  the  negative  thrust  limits  iTuposed  on  the  hookup 
simulations  were  much  greater—  -5000  or  -10000  lbs,  depending  on  the 
thrust  range  under  Investigation.  This  large  Imbalance  of  available 
forces  caused  no  problems  as  long  as  the  hookup  simulations  were  run 
under  no-gust  conditions.  As  a matter  of  fact,  the  hookup  simulations 
seemed  quite  realistic — the  pilot  model  put  in  the  full  500  lbs  of 
forward  thrust  at  the  precontact  position;  the  B'<S2  accelerated  slowly, 
attaining  a forward  velocity  In  the  neighborhood  of  1.5  to  1.8  ft/sec 
as  it  neared  the  refueling  envelope,  at  which  time  the  pilot  model  put 
In  a large  amount  of  negative  thrust  to  drive  the  closure  rate  to  zero 
as  the  aircraft  arrived  at  the  contact  position. 

When  wind  gusts  were  Introduced  into  the  B-52  equations  of  motion 
the  small  forward  thrust  limit  of  500  lbs  was  insufficient  to  allow 
the  pilot  model  to  adequately  correct  for  backward  movements  caused  by 
the  gust-induced  perturbations.  As  a result,  the  simulated  B-52  often 
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was  caused  to  move  backwards  to  the  limits  of  the  analog  scaling  of  the 
X-axls,  which  was  200  ft.  The  time  required  to  effect  a hookup  In 
these  Instances  was  usually  quite  long.  The  assumption  was  made  that 
the  pilot  would  not  maintain  the  low  positive  thrust  setting  (500  lbs) 
while  his  aircraft  drifted  backwards.  In  an  attempt  to  make  the  hookup 
simulation  more  realistic  In  the  wind  gust  environment,  a type  of 
graduated  forward  throttle  was  developed. 

By  means  of  analog  patching,  the  forward  throttle  limit  was  al- 
lowed to  Increase  above  500  lbs  if  a negative  velocity  along  the  X-axis 
was  experienced.  The  limit  increased  In  proportion  to  negative  u up  to 
Che  maxlmums  allowed  for  the  boom  track  and  Inadvertant  disconnect  pro- 
cedures (i.e.,  2500  or  5000  lbs).  The  thrust  limit  Imposed  by  this 
graduated  throttle  Is  given  by  the  schedule 


Forward  thrust  (lbs)  = 500  u ^ 0 

- 500 (1-u)  u < 0 


(57) 


Monte  Carlo  Simulation  of  the  Air  Refueling  Task 

Each  part  of  the  refueling  mission — hookup,  boom  track,  and  inad- 
vertanC  disconnect — was  simulated  using  the  ranges  of  variables  given 
In- Table  VIII.  Because  the  B-52  equations  of  motion  were  subjected 
to  gust  disturbances  of  a random  nature,  a series  of  500  refuelings 
were  simulated  for  each  combination  of  variables,  and  statistics  were 
compiled  on  the  duration  of  each  of  the  three  components  of  the  air 
refueling.  Each  of  these  simulations  was  timed  fot  2000  seconds  unless 
sooner  terminated  by  the  error  detection  and  logic  circuitry  on  the 
analog  computer. 
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For  the  boom  track  procedure  this  early  termination  occurred  at 
the  time  of  first  boom  disconnect  (l.e.,  when  the  B-52  boom  receptacle 
first  went  outside  the  limits  of  the  refueling  envelope). 

For  the  hookup  and  inadvertant  disconnect  procedures  the  runs 
were  terminated  in  less  than  2000  seconds  only  when  the  B-52  boom 
receptacle  had  remained  for  fifteen  seconds  within  the  hookup  and  pre- 
contact envelope,  respectively. 

All  of  the  simulations  were  run  on  the  analog  computer  at  a rate 
of  1000  times  real  time.  Thus  only  2 seconds  of  real  time  were  required 
to  simulate  2000  seconds  of  flight.  For  this  purpose  the  frequency 
bandwidth  of  the  random  noise  Input  to  the  gust  filter  also  was  In- 
creased by  a factor  of  1000,  to  1.5  KHz. 

At  the  termination  of  each  run  the  duration  in  seconds  for  that 
run  was  collected  and  stored  by  the  PDP-12  digital  computer.  After 
the  time  statistics  were  collected  for  each  series  of  500  runs,  the  f 

mean  time  for  that  series  was  calculated. 

Using  the  mean  times  for  hookup,  boom  track,  and  inadvertant  dls- 

i 

connect  the  mean  total  refueling  time  was  calculated  as  shown  below.  { 

RT  - Integer  x (HU  + ID)  + HU  + 1200  (58) 

i 

where  RT  * mean  refueling  time 

BT  ■ mean  tlme-on-boom  (boom  track  time)  j 

HU  ■ mean  hookup  time 

i 

ID  « mean  Inadvertant  disconnect  time 
1200  ■ total  tlme-on-boom  required  to  accomplish  the  refueling 

Integer  ■ Integer  value  of  which  represents  the 

number  of  inadvertant  disconnects  which  occur  during  the  refueling. 
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The  Initial  hookup  time  Is  accounted  for  In  the  total  time,  but  the 
final  disconnect  Is  not  Included  because  the  B-52  need  not  return  to 
the  pre-contact  position  and  stabilize  after  the  refueling  Is  completed. 
To  offer  a comparison  of  the  results  of  these  simulations  with 
the  results  had  the  downwash  from  the  KC-135  been  neglected,  a group 
of  runs  at  Tl  = 5 were  made  with  the  downwash  effects  excluded. 


V.  Results  and  Discussion 

Overview 

The  data  collected  In  the  Monte  Carlo  simulations  of  the  air 
refueling  task  are  presented  In  this  chapter.  The  data  are  In  the 
form  of  plots  of  the  trends  of  mean  refueling  time  as  functions  of 
rms  gust  intensity,  a„.  The  other  three  simulation  variables — thrust 
authority,  Tl,  and  radiation  dose — are  presented  as  parameters  on 
these  plots.  The  following  four  sections  will  discuss  the  variation 
of  mean  refueling  time  with  each  of  the  simulation  variables.  Finally, 
a brief  discussion  will  be  presented  on  the  results  of  the  simulation 
runs  with  no  tanker  downwash.  Appendix  H includes  supplementary 
plots  of  the  components  of  mean  refueling  time — hookup  time,  tlme-on- 
boom,  and  Inadvertant  disconnect  time — as  functions  of  gust  Intensity. 

Variation  of  Mean  Refueling  Time  with  Gust  Intensity 

Figure  7 shows  a plot  of  the  trend  of  mean  refueling  time  versus 
gust  intensity  for  the  non- irradiated  pilot.  Inspection  of  the  curves 
for  Thrust  Authority  A (narrow  range  of  thrust  authority)  reveals  that 
Che  mean  refueling  time  Increases  slowly  for  gusts  of  3 ft/sec  rms 
or. less.  Above  3 ft/sec  rms  the  refueling  time  increases  sharply, 
much  like  an  exponential  function  of  gust  intensity.  Similar  varia- 
tions are  exhibited  by  the  plots  for  500  and  1000  rads.  Figures  8 and 
9,  respectively.  The  mean  refueling  times  at  Thrust  Authority  B (wide 
range  of  thrust  authority)  Indicate  a sharp  rise  only  when  o„  exceeds 
5 ft/sec  rms. 


1 


II 


44 


0 1 2 3 4 5 6 


a^»  ft/sec  rms 


Figure  0.  Mean  Refueling  Time  vs.  Gust  Intensity  for  5OO  Rads, 
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a^,  ft/sec  rms 


Figure  9»  Mean  Refueling  Time  vs.  Gust  Intensity  for  1000  Rads. 


Variation  of  Mean  Rofuellnr.  Time  with  Thrust  Authority 

Figures  7,  8,  and  9 show  that  the  thrust  authority  given  the 
pilot  model  has  a pronounced  effect  on  the  mean  refueling  time.  As 
thrust  authority  Is  Increased  from  level  A to  level  B,  as  defined  In 
Table  VIII,  the  pilot  model  Is  able  to  perform  the  refueling  In  less 
time  at  the  higher  gust  Intensities. 

Variation  of  Mean  Refueling  Time  with  T^ 

Values  of  T^  are  shown  on  each  of  the  curves  on  Figures  7,  8,  and 
9.  Figure  7 shows  that  the  optimum  Tl  Is  5 seconds  for  the  non- 
Irradlated  pilot  model  operating  under  the  restrictions  of  Thrust 
Authority  A.  For  Thrust  Authority  B the  optimum  T^  appears  to  be  A 
seconds.  Figure  8 shows  that  Tl  4 sec  Is  the  optimum  pilot  lead 
for  the  case  of  a 500  rad  radiation  dose.  This  observation  holds  for 
both  levels  of  thrust  authority.  Figure  9 indicates  that  the  optimum 
Tl  for  1000  rads  also  Is  4 sec;  however.  In  the  case  of  Thrust  Author- 
ity A,  the  gap  between  4 sec  and  3 sec  appears  to  be  closing  somewhat. 
If  that  trend  continues,  3 sec  may  well  be  the  optimum  value  of  pilot 
lead  at  radiation  doses  greater  than  1000  rads. 

Variation  of  Mean  Refueling  Time  with  Radiation  Dose 

Perhaps  the  most  interesting  effect  on  this  model  of  the  air 
refueling  task  Is  that  of  pilot  irradiation  as  simulated  by  gain  decre- 
ments. Examination  of  Figures  7,  8,  and  9 shows  that  there  is  rela- 
tively little  difference  in  mean  refueling  time  as  the  simulated 
radiation  dose  is  increased.  Two  other  trends  assume  greater  im- 
portance here.  One  of  these  trends-*- the  shifting  of  the  optimum  Tl 
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to  a lower  value  as  radiation  dose  increases — was  discussed  in  the 
previous  section.  The  other  trend  is  clearly  shown  only  by  the  plot 
for  1000  rads.  Figure  9.  The  curves  for  Tl's  equal  to  5 and  6 sec 
were  not  run  because  with  Tl  equal  to  6 sec  the  pilot-aircraft  system 
was  unstable  in  the  elevator  control  loop,  and  the  refueling  time  data 
collected  when  Tj^  was  set  equal  to  5 sec  was  erratic  to  the  point  of 
being  unuseable.  The  values  of  Kp^  used  in  these  Instances  were 
.003281  and  .003482  for  T^  equal  to  6 and  5,  respectively.  These  pilot 
gains  were  derived  by  multiplying  the  non-irradlated  pilot  gains  given 
in  Table  V by  a gain  factor  of  .83  (Table  VI).  Table  V also  shows 
that  the  minimum  stable  Kp  is  .003404  for  Tj^  = 6 and  .003403  for 
Tl  = 5 sec.  Clearly,  the  value  of  1C,  used  in  conjunction  with 
■ 6 was  in  the  unstable  region  while  the  value  of  used  with 
Tl  “ 5 was  above  but  very  close  to  the  minimum  stable  value.  This 
marginally  stable  aspect  of  the  pilot-aircraft  system  having  Tl  equal 
to  5 sec  accounts  for  its  erratic  behavior  in  the  simulations. 


Results  of  No-Downwash  Simulations 

Figure  10  depicts  the  trend  of  mean  refueling  time  with  Increasing 
gust  Intensity  for  the  model  of  the  non-lrradlated  pilot  with  Tl  5 
sec.  The  data  for  both  levels  of  thrust  authority  show  marked  In- 
creases In  refueling  time  over  the  corresponding  curves  with  downwash 
that  are  presented  In  Figure  7.  The  conclusion  may  have  been  drawn 
from  the  comparisons  of  Tables  III  and  IV  that  the  downwash  effects 
caused  the  refueling  task  to  be  more  difficult  due  to  decreasing  the 
damping  and  Increasing  the  frequency  of  the  short  period  mode.  These 
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results  Indicate  that  the  air  refueling  task  Is  actually  made  less 
difficult  by  the  Introduction  of  the  downwash.  This  phenomenon  was 
not  Investigated  further  than  to  observe  that  the  rise  time  of  the 
system  with  downwash  was  shorter  than  that  of  the  system  without 
downwash.  In  any  case,  the  no-downwash  situation  is  not  encountered 
in  practice. 
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VI.  Conclusions  and  Recommendat Ions 


The  preceding  discussion  showed  that  the  trends  displayed  by  the 
mean  refueling  time  as  gust  Intensity  increased  depended  heavily  on 
the  limits  of  the  pilot  thrust  authority  imposed  on  the  simulation. 

It  is  not  clear  at  this  point  as  to  how  much  thrust  authority  can 
realistically  be  allowed  the  pilot  model  in  the  simulation  of  the 
air  refueling  task.  It  is  recommended  that  this  aspect  of  the  problem 
be  Investigated  In  future  work,  to  Include  simulator  studies  involving 
human  subjects  and  the  collection  of  field  data  during  actual  air  re- 
fueling missions. 

Simulation  of  the  air  refueling  task  at  two  levels  of  thrust 
authority  provided  some  insight  into  the  variation  of  Tl  with  diffi- 
culty of  task.  The  data  seemed  to  indicate  that  the  optimum  value 
of  Tl  tended  to  be  lower  for  the  runs  with  Thrust  Authority  B,  the 
large  range  of  thrust  authority.  Although  this  trend  was  subtle,  one 
might  assume  that  the  optimum  Tl  decreases  as  the  difficulty  of  the 
task  decreases. 

The  conclusion  that  seems  to  be  indicated  by  the  results  of  the 
simulated  pilot  irradiation  by  means  of  gain  decrements  is:  the  pilot 

must  decrease  his  Tl  to  maintain  stability  as  his  radiation  dose 
increases.  If  the  gain  decrement  causes  Kq  to  decrease  below  the 
minimum  stable  value,  there  is  no  recourse  to  the  pilot  except  to 
lower  his  Tl  to  assure  system  stability. 

Much  work  needs  to  be  done  to  refine  the  model  of  the  air  refuel- 
ing task.  An  important  step  would  be  the  inclusion  of  the  lateral 
equations  of  motion.  This  step  would  require  the  development  of 
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APPENDIX  A 


Linearization  of  Induced  Effects  Data  by  the  Method  of  Least  Squares 

Table  IX  presents  the  Induced  effects  of  the  KC* 135  flow  field 
(downwash)  on  the  lift  and  moment  coefficients  of  the  B-52  in  the 
refueling  situation. 

Table  IX 

Induced  Effects  Data 


Inches 

ACli 

-39 

-.080 

.072 

11 

-.077 

.067 

61 

-.072 

.063 

111 

-.066 

.058 

161 

-.061 

.054 

211 

-.057 

.050 

261 

-.052 

.046 

311 

-.048 

.042 

(Ref  5:  2-79,4-95) 


Method  of  Least  Squares: 

To  fit  the  linear  equation  y = aQ  + aj^x  to  a set  of  data  (Xj.Yj^), 
1 * l,n  It  Is  necessary  to  solve  for  the  constants  ag  and  aj^  by 
means  of  the  simultaneous  equations: 


EYi  - ag  n + ai  EX^^  (59) 

EXjYi  - ag  EXjl  + ai  EX^^ 


(60) 


To  find  as  a function  of  z:  = Zj.  (61) 

Yi  * (ACLj^)i  (62) 

n ■ 8 (63) 

The  following  simultaneous  equations  are  formulated: 

Sag  + 1088ai  = -0.513  (64) 

1088ao  + 252968ai  » -59.793  (65) 

Solving  yields  Hq  • 0.000095  (66) 

ai  » -.077045  (67) 

Therefore,  the  linear  equation  for  is 

- -.077045  + 0.000095Z  (z  in  inches)  (68) 

or  * -0.077045  + 0.00114z  (z  in  feet)  (69) 

Similarly  for  the  simultaneous  equations  are 

8ao  + 1088ai  » 0.452  (70) 

1088ao  + 252968ai  = 52.522  (71) 

which  yield  Sq  « 0,0680835  (72) 

ai  - -0.0000852  (73) 

The  linear  equation  for  is 

* 0.0680835  - 0.0000852z  (z  in  inches)  (74) 

or  * 0.0680835  - 0.0010224z  (z  in  feet)  (75) 
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APPENDIX  B 


Calculation  of  Steady  State  Angle  of  Attack  In  Dovnwash 

In  the  case  of  steady  state  flight  with  no  downwash,  the  follow- 
Ing  equations  are  applicable: 


E °stab 

<1.1  - ° ■ q.  + q.  ->1  q..  h*  «.t.bi  <”> 

° “ E * *stab 

If  the  aircraft  subsequently  Is  trimmed  for  steady  state  flight  In  a 
downwash  flow  field  the  steady  state  equations  are 


Cl^  + + Ojjy)  + + C 


Lfis^ab 

+ ACt,  (78) 


0 - Ca^  + Cn^  (Oj^  + (^stabi  + ^stabi,„^ 

+ ACa^  (79) 

frtiere  Is  the  horizontal  stabilizer  input  which  produces  the 

additional  angle  of  attack,  Oj^y,  required  to  trim  out  the  downwash- 
Induced  changes  In  the  lift  and  moment  coefficients,  sod  Cq^, 
respectively. 

If  the  weight  of  the  aircraft  does  not  change,  then  must 
remain  the  same  before  and  after  entering  the  downwash  field.  The 
moment  coefficient,  Caj^>  must  also  remain  zero  if  the  aircraft  is  in 
equilibrium  level  flight.  Therefore 


CLa  «DW  + «8tabDW  ° 

Cma  “dW  + ‘^Astab  ^stabow  + ° 
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For  trlnuned  flight  at  the  center  of  the  refueling  envelope  (z  =•  0) 
these  equations  are 


( .0891)  ajjj,  + (.01079)  «stabD„  " .077  = 0 (82) 

(-.0151)  ojjw  - (.0327)  SstabDW  ^®3) 

Solving  yields 

Ojjj,  - 0.65°  (84) 

^stab^y  * 1.78°  (85) 

The  steady  state  angle  of  attack  required  to  maintain  the  B-52 
in  equilibrium  flight  at  the  center  of  the  refueling  envelope  is 

oj  (in  downwash)  = 0.03°  + 0.65°  * 0.68°  (86) 


59 


APPENDIX  C 


Details  of  the  Derivation  of  the  Controlled  Elements 

The  required  aircraft  transfer  functions  are  derived  from  Eq  (44) 
by  means  of  Cramer’s  Rule.  The  gust  terms  are  ignored  in  these 
derivations.  A and  B continue  to  be  used  to  represent  the  coefficients 
of  2,  the  dovmwash  variable.  The  transfer  functions  which  are  needed 
for  the  derivation  of  the  controlled  elements  are 


-8.143  - 5.3384  + 166.705s  - ,071184 

tTE: 


e -.627  --239  +(-.0025167+.627A-8.143B)S+(.003A-.036B) 


. 3 -i 

«E  «E 


C.E. 


(87) 


(88) 


(89) 


^ _ .0000715  (s ^+1.417  33+(2.291-A)s2+(-1.005A+4.6B)s+(-1.9A+290.6B)  } 
«T  C.E. 

where  C.E.  is  the  general  characteristic  equation  as  given  by  Eq  (45). 
Eqs  (87)  and  (88)  above,  are  then  used  to  calculate  by  the  method 


of  Eq  (12).  That  calculation  yields 


^ _ 25.7  s3  + 7.6  s2  ■Kl66.8-34A+440B)s-K-.0712-.15A+2B) 
6g  C.E. 


(91) 


By  block  diagram  reduction  of  the  B-52  pitch  SAS  control  loop  shown 


in  Figure  2,  the  transfer  function  was  found  to  be 

*pe 


.9315  s(s+.15)(s-b.25)(^/*E) 


®PE 


(92) 
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where  A - s^+2.45s^+(3.14625-A)s5+(. 988-2. 0A+12.8B)sV(. 099-2. 3A+292.6B)s3 

+(.0037-.8A+116.4B)s2+(.0003-,074A+10.9B)s+(-.00032A-.000057B) 

(93) 

Finally,  the  elevator  controlled  element  Is  derived  by  substituting 

Eqs  (89),  (91),  and  (92)  into  Eq  (48),  yielding  Yj.  = k — 

^ °PE 

_ .93(s+.15)(s+.25){25.7s3+7.6s2+(166.8-34A+440B)s+(-.0712-.15A+2B)} 

A 

(94) 

The  throttle  controlled  element  is  derived  by  substituting  Eqs  (25) 
and  (90)  into  Eq  (49) . 

The  controlled  elements,  with  and  without  downwash,  are  presented 


below  in  factored  form. 


Without  Downwash:  A = B = 0 


y ^ 23. 93(s+. 15) (s+. 25) (s+. 1475±1 2 . 54) (s-. 0004)  - . 

No  DW  " s(s+1.03+jl.ll)(s+.0017+j.0593)(s+.2216)(s+.1673) 

Y ) m .0000715  (96) 

No  DW  (. 25S+1)  (s+. 0017+j . 0593)  ^ 


With  Downwash: 


-.06471,  B - -.002269 


y'  ) 

DW 


^CT^DW 


23 . 93  (s-t-.  1474+1 2 . 553)  (s- . 0004) 

(s+. 9066±jl. 158) (8+.631)(s-. 3891) (s-. 00105) 

. 0000715 (s+. 54945) (s-. 40764) 
s(. 258+1) (S+.5484) (8-. 402) (8-. 00106) 


Bode  magnitude  and  phaae  angle  plot8  of  the  controlled  eleraenta 
without  downwaah  are  given  in  Figures  C-1  and  C-2.  Bode  plots  of  the 
controlled  elements  with  downwash  are  shown  in  Figures  C-3  and  C-4. 
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Figure  C-3.  Bode  Plot  c 
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APPENDIX  D 


Gust  Model 


Assume  the  Dryden  spectral  form 


— , ^ /3  Lw 

^ ^ Ul 


wg(s)  _ / Ul 

“ C(s)  ■ iruY  (1  + ^ s)2 
According  to  the  Dryden  scales  for  clear  air  turbulence, 


Then, 


Lw  = 1750  ft 


”g(s)  0 /T750 

5(s)  “ it(706) 


, . /3  1750 

(1  s)2 


„ ( OORON  1 + 4.29  s . 

®w  l-oo8i)  ^ ^ 2.48  s)' 


(100) 


where  a„  is  the  rms  gust  Intensity. 


(Ref  7:  444-459) 


Calculation  of  Gain  Needed  to  Compensate  for  Filter  Loss 
Assume  the  following  system; 


If  P(ui)  and  Q(u»)  are  power  density  functions  then  a filter  G(jw)  will 
have  the  effect: 


Q(oj)  - iG(ju)  p P(u)) 


66 


Therefore,  the  total  power  In  Q(u))  across  all  frequencies  Is; 


Q = /”Q(a))  do)  = /”|G{jo))|^  P(o))  dw 


(102) 


or 


q = /‘*|G(j0))12  P(o,)  do) 


(103) 


where  = 20  db  rolloff  point  of  the  filter  or  the  20  db  rolloff 

point  of  P((i)),  whichever  is  less.  In  the  case  of  the  gust  filter. 


\ 1 + a iu 

G(ju,)  Ow 


(104) 


where  a 


Ul 


Assume  uniform  power  density  to  1.05  fg,  where  fg  is  the  cutoff 
frequency  of  the  noise  generator. 


Then 


P(o>) 


2nfg(1.05)  2.1irfg  Old 


(105) 


Now  we  have: 


|G(ju))l 


2 . “ a 2 


1 + 3 0)^ 

(1  + g2  0)2)^ 


and 


(106) 


o 0, 


_ **  "'w  P 1 + 3 g~  U)~ 

^ ” IT  ud  0 (1  + g2  (j2)2 


1 + 3 


du 


(107) 


Q - - J-  a ^ (0) 
^ IT  Ud  W ' 


(108) 


where 


n ,"d  1 + 3 , 

i (1  + 0)^)2 


or 


^ dw  . ^ 2 Mi)A 

TrT^2^2  + 3o  — + cx2  ^2)2 


da> 


+ 3 o2  £22 


Integrating  yields 


“d  1 

* 2(1  -Ta^  a.  Z)  + 2^  ‘^’d) 

d 

1 ( -^d  1 

^2  - ^ i-+  C.-2  ^,2  + a arctan  (a 


For  a cutoff  frequency,  fg,  equal  to  1.5  Hz, 


(oj  ■ 2.1irfQ  = 9.896  rad/sec 


Substituting  the  numerical  values  for  Lw  and  U]^  yields 


a - 1750/706  - 2.48 


Then  £2]^  and  £22  are  found  to  be 


- 0.3166 
(22  - 0.0488 

Substituting  the  values  of  (2^^  and  (^2  Into  Eqs  (110)  and  (108) 

Q - o„2  {0.3166  + 3(2. 48)2(. 0488)  } 

- .09708  o„2  p 

i/q"-  /.09708  /T 


(109) 

(110) 

(111) 

(112) 

(113) 

(114) 

(115) 

(116) 

yields 

(117) 

(118) 
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Therefore,  the  gain  needed  to  make  up  for  the  filter  loss  is 


1 ^ 3.209 

/.09708 


Kp  is  applied  to  the  gust  model  as  an  additional  factor  in  the 
so  that  Eq  (100)  becomes 


^e(s) 

5(s) 


2.85 


4.29  s + 1 
w (2.48  s + 1)2 


which  is  the  gust  model  (Ref  13). 

A bode  plot  of  Eq  (120)  is  shown  in  Figure  D-1. 
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APPENDIX  E 

Bode  and  Polar  Plots  of  the  Open-Loop  Pilot-Aircraft 
Transfer  Functions  Without  Downwash 


Bode  and  polar  plots  of  Yp^  and  Yp^  are 

shown  for  values  of  Tl  starting  with  Tl  = 8 and  decreasing  to  Tl  = 1. 
These  transfer  functions  are  given  by  Eqs  (52)  and  (53).  Pilot  gains 
used  in  calculating  these  plots  ate  given  in  Table  V. 
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Figure  E-6,  Bode  and  Polair  Plots  of  Y Y L.  T,  - 6, 
^ Pj  c^'No  DW  L 
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APPENDIX  F 


Bode  and  Polar  Plots  of  the  Open-Loop  Pilot-Aircraft 
Transfer  Functions  with  Downwash 

The  figures  contained  herein  present  the  Bode  and  polar  plots 

of  Yn  Y-  ) and  Yn„  ^Co^^nu  Tt  Is  varied  from  8 to  1.  The  transfer 
E DW  ^T  ‘^T  DW 

functions  are  given  by  Eqs  (54)  and  (55).  The  values  of  K„  used  in 
Eq  (54)  are  listed  in  Table  V.  The  plots  of  Yp^  ^c^^DW  were 
calculated  by  setting  Kp^  equal  to  1/. 06857  = 14.58  in  Eq  (55)  in 
order  to  show  more  detail  in  the  -1  + jO  regions  of  the  polar  plots. 
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APPENDIX  G 


Analog  Patching  Diagram 

The  scaled  analog  patching  diagram  Is  shown  In  Figure  G-1. 
Table  X lists  the  potentiometer  settings. 


Table  X 

Analog  Potentiometer  Settings 


Pot 

Setting 

Pot 

Setting 

Pot 

Setting 

0 

.0045 

32 

.2500 

74 

.1000 

1 

.0500 

33 

.2500 

75 

.3447 

2 

.0091 

34 

.0141 

76 

.1990 

3 

.3220 

35 

.2500 

90 

.3333 

4 

.4117 

36 

.1863 

91 

.3333 

6 

.0715 

37 

.0108 

92 

Tl/10 

.2548 

38 

.1082 

94 

.3779  hookup 

7 

.3367  boom  track 

8 

.4117 

40 

.1082 

.8500  disconnect 

9 

.0701 

41 

.0407 

96 

.3333 

10 

.0906 

42 

.0314 

97 

.2000 

11 

.0034 

43 

. 1863 

98 

Kp^/ 10000 

12 

.7060 

46 

.0540 

99 

.4663  hookup 
.3640  boom  track 

13 

.2734 

50 

.0540 

.2679  disconnect 

14 

.5400 

51 

.2734 

101 

.4000 

15 

f 

.7060 

53 

.2548 

102 

.4000 

16 

.0228 

60 

.3333 

103 

.1000 

18 

.0647 

61 

.3333 

110 

.5000 

19 

.3209 

62 

T^/IO 

111 

.1414 

20 

.2331 

64 

.4179  hookup 
.4592  boom  track 

Q2 

.2500 

21 

.8065 

.9500  disconnect 

Q7 

.1000 

22 

.6975 

66 

.3333 

Q9 

.2500 

23 

.6196 

67 

.2000 

Q17 

.1000 

30 

.1500 

68 

o 

o 

n 

Q24 

o„/10 

31 

.1500 

69 

.7002  hookup 
<8391  boom  track 

.4663  disconnect 


APPENDIX  H 

Plots  of  Components  of  Mean  Refueling  Time 

The  components  of  mean  refueling  time— hookup  time,  time-on- 
boom,  and  Inadvertant  disconnect  time — are  presented  here  as  func- 
tions of  gust  Intensity  with  T^,  radiation  dose,  and  thrust  authority 
as  additional  parameters.  The  plots  of  Inadvertant  disconnect  times 
depict  sketchy  Information  due  to  the  fact  that  the  Inadvertant  dis- 
connect procedure  usually  was  simulated  only  when  required  for  Eq  (58). 


- 


Mean  Hookup  Time 
(Trend) 


Figure  H-2.  Mean  Hookup  Time  vs.  Gust  Intensity  for  0 Hauls 
Thrust  Authority  B. 


Mean  Time-on-Boom 
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Flgxire  H-3.  Mean  Time-on-Boom  vs.  Gust  Intensity  for  0 Rads 
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4^ 


Figure  H-6.  Mean  Hookup  Time  vs.  Gust  Intensity  for  5OO  Rads, 
Thrust  Authority  B. 


Mean  Time-on-Boom 


Note:  values  of 

li 

shown  on  curves 


T - I 1 1 1 1 1 

0 1 2 3 4 5 6 

cr^f  ft/sec  rms 


Figure  H-8.  Mean  Inadvertant  Disconnect  Time  vs.  Gust  Intensity 
for  500  Rods. 
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Figure  H-9.  Mean  Hookup  Time  vs.  Gust  Intensity  for  1000  Rads, 
Thrust  Authority  A. 
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Figure  H-10,  Mean  Hookup  Time  vs.  Gust  Intensity  for  1000  Rauls • 
Thrust  Authority  B. 
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Mean  Tlue-on-Boom 


^ Figure  H-I3.  Mean  Hookup  Time  vs.  Gust  Intensity  for  0 Bads» 

No  Downwush. 


125 


Note*  5 


Thrust 
Authority  A 


Thrust 
Authority  B 


a^t  ft/ sec  rms 


Figure  H~14.  Mean  Tlme-on-Boom  vs.  Gust  Intensity  for  0 RauiSf 
No  Downwash. 
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Figure  H-15t  Mean  Inadvertant  Disconnect  Time  vs.  Gust  Intensity 
for  0 RadS(  No  Downwash. 
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The  controlled  elements  with  and  without  downwash  are  derived  for  the 
refueling  task.  Comparison  of  these  transfer  functions  shows  that  the 
downwash  causes  Instability.  Stability  analyses  are  conducted  on  the 
pilot-aircraft  open-loop  trMsfer  functions  over  a range  of  pilot  rate 
to  displacement  weighting,  a , and  pilot  gains  are  found  which  will 
stabilize  the  system  for  eactr  value  of  . '^T 

Monte  Carlo  simulations  of  the~~fef uel^g  task  are  run  at  varying 
combinations  of  , rms  gust  intensity,  pilot  thrust  authority,  and 
pilot  gain  decrement  due  to  Irradiation.  Trends  of  mean  refueling  time 
are  predicted.*^  Using  mean  refueling  time  as  a basis  for  comparison, 
optimum  values*^!  T are  found  for  the  Irradiated  and  the  non-irradiated 
pilot. 
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